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Abstract — We theoretically study the magnetoresistance of a 
CPP-GMR system with current confined paths (CCP) in the 
framework of Valet-Fert theory. The continuity equations for 
charge and spin currents are numerically solved with the three- 
dimensional CCP geometry by use of finite element method. It is 
confirmed that the MR ratio is enhanced by the CCP structure, 
which is consistent with the experimental results. Moreover, we 
find that there exists a certain contact width which maximize the 
MR ratio. We show that the contact width which maximize the 
MR ratio is well described by the effective resistance matching. 

Index Terms — CPP-GMR, current-confined-path, spin accu- 
mulation, nano-oxide-layer 

CURRENT-perpendicular-to-plane giant magnetoresis- 
tance (CPP-GMR) has attracted much attention for its 
potential application as a read sensor for high-density magnetic 
recording IT], ||2j. In order to realize a high-density magnetic 
recording, we need MR devices with high MR ratio and low 
resistance area product (RA). Although the RA value of a 
CPP-GMR system is much smaller than that of a tunneling 
magnetoresistance (TMR) system, the MR ratio of a conven- 
tional CPP-GMR system still remains a small value of a few 
%. Much effort has been devoted to increasing the MR ratio 
of the CPP-GMR system. 

Recently, Fukuzawa et al. reported that they achieved the 
MR ratio of 10.2 % by CPP-GMR spin-valve with a current- 
confined-path (CCP) structure made of a nano-oxide-layer 
(NOL) with a lot of small metallic channels 13J, ||4J. They 
showed that the MR ratio increases with increasing the RA 
value, which implies that the MR ratio is enhanced for the 
system with narrow metallic channels. The similar enhance- 
ment of the MR ratio due to the CCP structure containing 
a domain wall was also reported by Fuke et al. 0. The 
enhancement of the MR ratio due to the CCP structure is 
theoretically explained by one of the author |6|. However, the 
analysis in |6J is based on the assumption that the thickness 
of the non-magnetic layer is zero and little attention is paid to 
the effect of the nonmagnetic spacer layer on the enhancement 
of the MR ratio. 

In 2000, Schmidt et al. pointed out that the conductivity 
mismatch is a basic obstacle for spin injection from a ferro- 
magnetic metal into a semiconductor |7|. Spin polarization of 
the injected current is strongly suppressed by the conductivity 
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mismatch between these materials. For conventional CPP- 
GMR spin-valves without CCP structures, the conductivity 
mismatch seems not to be important because both of the 
ferromagnetic electrodes and the nonmagnetic spacer layer 
are made of metals. However, for the CCP-CPP-GMR system 
the effective resistance of the contact region increases with 
decreasing the radius of the contact, which implies that the 
conductivity mismatch plays an important role in the CCP- 
CPP-GMR system. Therefore, it is intriguing to ask how the 
mismatch of the effective resistance affects the MR ratio of 
the CCP-CPP-GMR system. 

In this paper, we analyze the dependence of the MR ratio 
of the CCP-CPP-GMR spin-valves on the contact radius and 
resistance area product (RA) in the framework of Valet-Fert 
theory [8|. The spin-dependent electro-chemical potentials 
are obtained by numerically solving continuity equations for 
charge and spin currents with the three-dimensional CCP 
geometry by use of the finite element method ||9l, ifTOl . We 
show that the CCP-CPP-GMR spin-valve can take a larger 
MR ratio than that of the conventional CPP-GMR spin-valve. 
However, the MR ratio is not a monotonic function of the 
contact radius but takes a maximum value at a certain value of 
the contact radius. We also show that the contact radius which 
maximizes the MR ratio can be explained by considering the 
matching of effective resistances. 

I. Valet-Fert Theory of CPP-GMR 

By starting with the Boltzmann equation. Valet and Fert 
constructed the macroscopic model of the CPP-GMR|8|, on 
which our calculation is based. Before showing our results, 
we shall give a brief introduction to the Valet-Fert theory of 
CPP-GMR. 

In the macroscopic model of the CPP-GMR, it is assumed 
that the conduction electrons with spins s =t, j can be charac- 
terized by the conductivity a-g, spin relaxation rate r^, density 
of states Ng and diffusion constant Dg. Since we consider 
a collinear alignment of the magnetizations we neglect the 
mixing of spin-up and spin-down bands. The current density 
for each spin band under an electric field E is given by 



asE - qDsVSus 



(1) 



where q is the charge of an electron and dug represents the 
accumulation of the electrons with spin s. By using Einstein's 
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relation — q^NsDg, we can rewrite ([TJ as 



(a) 



(b) F 



= --[-qE 

q 



(2) 



Introducing the deviation of the spin-dependent chemical 
potential from the equilibrium state 5 jig, the accumulation 
Sus can be expressed as Sus — NsSfis- Since the electric 
field is given by the gradient of the electrostatic potential cj) 
as E — ~'V(f), it is convenient to define the spin-dependent 



electro-chemical potential as fis = Sjis 
density ([T]) can be written as 

3s = V^s 



Hence the current 



(3) 



In the ferromagnetic materials, the spin-polarization of the 
conductivity is represented by the parameter (3 which is 
defined as (3 — ((T| — (Tx)/((T| + cj.)- For example, the spin- 
polarization parameter (3 = 0.36 ^ 0.5 for Co ifTTI 

The equations that determine electro-chemical potentials are 
derived by the continuity equations for charge and spin. The 
continuity equation for charge is given by 



V-(jT+Ji) =0, 



(4) 



which means that there is no source for electron charge. The 
continuity equation for spin is given by 



V-(jT-Ji)-9( — Mi), (5) 

which means that the spin relaxation plays a role of spin 
source. Substituting (O into ^ and ^ we obtain the fol- 
lowing diffusion equations for electro-chemical potentials 



cr-fV /i| 

V' (mt - 



f dxVVx = 0, 

Mi) ^(mt - 



MiJ 



(6) 

(7) 



where A is the characteristic length of the spin diffusion called 
"spin diffusion length". The spin diffusion length is defined as 
A — VtD, where t is the averaged spin relaxation time de- 
fined as — (r|~^+Tj~^)/2 and D is the averaged diffusion 

constant defined as (^T^j^^ + ^i^f ^)/(^T + ^i)- 

By solving ^ and (|7]i with proper boundary conditions, 
we obtain the electrochemical potential and therefore the total 
resistance of the system. The MR ratio is defined by the total 
resistance for the parallel alignment of the magnetizations Rp, 
and that for the anti-parallel alignment R^p as MR = (-Rap ^ 

II. CPP-GMR WITHOUT CCP 

Let us first consider the MR ratio of the conventional CPP- 
GMR spin-valves shown in Fig. [T] (a). The nonmagnetic (N) 
spacer layer is sandwiched between two ferromagnetic (F) 
electrodes. The magnetization of the top electrode is assumed 
to be fixed at a certain direction and the magnetization of the 
bottom electrode is aligned to be parallel or anti-parallel to 
that of the top electrode. 

We assume the translational invariance in the x— and 
y— directions and the current is flowing perpendicular to 
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Fig. 1. (a) A conventional CPP-GMR is scliematically shown. The 
nonmagnetic (N) spacer layer is sandwiched between two ferromagnetic (F) 
electrodes. The cun'ent is flowing perpendicular to the plane. The magneti- 
zation of the top electrode is assumed to be fixed at a certain direction and 
the magnetization of the bottom electrode is aligned to be parallel or anti- 
parallel to that of the top electrode, (b) Equivalent circuit diagram representing 
the total resistance of the conventional CPP-GMR spin-valve is shown. The 
resistance due to spin accumulation AR can be expressed as a parallel circuit 
of the effective resistances, (c) MR ratio is plotted against the resistivity of 
the non-magnetic layer pjv. (d) MR ratio is plotted against the ratio of the 
effective resistances r^/fjv- 



the plane, i.e., in the z— direction. Hence the equations to 
be solved reduce to a couple of one-dimensional diffusion 
equations: 



^ig^lMi 



0, 



1 

^(mt - Mi) = i^(mt - Mi)- 

The general solutions of (H) and ^ are obtained as 

1 



Mi (2) 



Cl + C2Z 



Cl + C2Z - 



1 

^i 



z/X 



-z/\ 



€46 



(8) 
(9) 

(10) 
(11) 



where c/s are arbitrary integration constants determined by 
the proper boundary conditions. 

We neglect the effect of interfacial scattering for simplicity. 
Then we adopt the boundary conditions such that both /i^ 
and is are continuous at the interface. After some algebra, 
we can show that the total resistance for parallel alignment of 
the magnetizations Rp, and that for the anti-parallel alignment 
Ra?, have the form 



-RrAP — -Ro + Ai?pAP, 



(12) 



where i?o = {^tppp + tNPN) /S, <i?(jv) is the thickness of 
the F(N) layer, S is the cross-section area of the system, and 
Pf{n) is the resistivity of the F(N) layer defined as 

1 1 

pp — ~f — : — r Miv — — -■ (13) 



F ~ " F 

' F(N) i^ '■^^ conductivity for electrons with spin s in 
the F (N) layer. Rq is the total resistance of the system 



Here ai 
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Fig. 2. (a) A CCP-CPP-GMR spin-valve is scliematically shown. Tlie 
nonmagnetic (N) contact region is sandwiclied between two ferromagnetic (F) 
electrodes, (b) The cross-sectional view of the CCP-CPP-GMR spin-valve is 
schematically shown. The curve of the contact region is modeled by the half 
ellipse, which is tangent to the ferromagnetic layer. In our calculation we keep 
the shape of the ellipse fixed and vary the size of the contact by changing the 
contact radius d. 



in the absence of the spin accumulation and independent of 
the magnetization configuration. The additional resistance AR 
originates from the voltage drop at the interfaces due to the 
spin accumulation, which has the form 



' N 



y 



(14) 



where rp and rjv are the effective resistances for F and N 
layers defined as 



5(1-/32)' 



S 



/p. 



2A 



N 



(15) 



The function /pap is given by fp{x) = tanh(2;) and fAp{x) = 
coth(a;). From ( fT4] i. AR is considered as a parallel circuit of 
the effective resistances rp and tn as shown in Fig. [r|(b). 

In Fig. [T] (c), we plot the MR ratio against the resistivity 
of the non-magnetic layer p^. We assume that the F layer is 
made of CoFe and the N layer is made of Cu. For numerical 
calculation, we use the following values: pp = 160ilnm, 
tp — 45 nm, — 2nm, Xp = 15 nm, \n ~ 500 nm, 
/3 = 0.7. The thickness of the F-layer tp is taken to be thick 
enough compared with the spin diffusion length Af in order 
to eliminate the spin accumulation at the bottom and the top 
electrodes. 

One can obviously see that the MR ratio takes its maximum 
value at a certain value of pjv and vanishes in both limits of 
Pjv — > oo and 0. In the limit of oo, the effective 

resistance of N-layer rjv is much larger than that of F-layer 
rp and the interface resistance AR equals to 20^rp for both 
P and AP cases, which means that the MR ratio goes to zero. 
Conversely, in the limit of pjv 0, rjv is much smaller than 
rp and the interface resistance AR goes to zero. Therefore 
the MR ratio vanishes also in this limit. 

In Fig. [T] (d), the MR ratio is plotted against the ratio of 
the effective resistance rp/f^, where fpf — pn^n is the 
effective resistance of N-layer for an F/N bilayer system. As 
shown in Fig.[T](d), the MR ratio is maximized if the effective 
resistances rp and fjv are balanced with each other, i.e., rp ~ 
fN- A similar idea is also applicable to the CCP-CPP-GMR 
system as will be discussed in the next section. 
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Fig. 3. MR ratios of CCP-CPP-GMR spin- valves are plotted against the 
contact radius d. The resistivity of Cu in the contact is varied from 17 Qnm 
to 300 Qrnn. The MR ratio takes a maximum value at a certain value of the 
contact radius which is indicated by a circle. 



III. CCP-CPP-GMR 

Next we move onto the CCP-CPP spin-valves, where the 
spacer layer is made of an NOL with a lot of nanometer-size 
metallic channels. The system we consider is schematically 
shown in Fig.|2](a). A narrow non-magnetic metallic contact is 
sandwiched between two ferromagnetic electrodes. The cross- 
sectional view of the system is shown in Fig.|2](b). The curve 
of the contact region is modeled by the half ellipse, which is 
tangent to the ferromagnetic layer The length of the major 
axis of the ellipse, which is the same as the thickness of the 
contact, is taken to be 2 nm. The length of the semi-minor 
axis of the ellipse is 0.5 nm. The F electrodes are assumed to 
be 40 nm x 40nm x 45 nm rectangles. 

In our calculation we keep the shape of the ellipse fixed and 
vary the size of the contact by changing the contact radius d. 
The electro-chemical potentials are determined by numerically 
solving ^ and (|7]l by use of the finite element method |]9], 
|10|. The system is divided into hexahedral elements and the 
total number of the elements is of the order of 10^. We assume 
that the F electrodes and N contact are made of CoFe and Cu, 
respectively. We use the same material parameters we used 
in the previous section for the conventional CPP-GMR spin- 
valve system. The resistivity of Cu is assumed to be pcu =17, 
50, 170 and 300 num. 

The obtained MR ratios of the CCP-CPP-GMR system 
are plotted against the contact radius d in Fig. [3] One can 
immediately see that the MR ratio increases with decreasing 
the resistivity of Cu, pcu- It has been considered so far that 
the MR ratio shows monotone increasing with decreasing the 
contact radius both theoretically |6| and experimentally [4). 
We note that there is a contact radius that maximize the MR 
ratio, which is indicated by a circle in Fig. [3] for given values 
of pcu and the MR ratio vanishes in the limit of d ^ 0. Fig. |3] 
also shows that the value of the contact radius that maximize 
the MR ratio decreases with decreasing the resistivity of Cu. 

Strictly speaking, the Boltzmann approach loses its validity 
in the limit d ^ 0, especially in the region d < Inm. 
In Ref. fT2^, the electron transport across a narrow metallic 
channel with the width of a single atom is studied based 
on the Landauer formalism, where it is concluded that MR 
ratio is diminished by the quantum interference. This does not 
contradict our results that MR goes to zero in the limit d ^ 0. 
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Fig. 4. (a) MR ratios of CCP-CPP-GMR spin- valves sliown in Fig. [3] are 
plotted against the resistance area product (RA). (b) MR ratios sliown in Fig. 
[5] ai'e plotted against the ratio of the effective resistances rp/fpf- In both 
panels, circles indicate the maximum value of the MR ratio for each value of 
PCii- 

In Fig. |4] (a), we plot the MR ratios against the resistance 
area product (RA). This behavior is quaUtatively consistent 
with the experimental results [31 ,0. The difference of the 
absolute scale of the MR ratio between our theory and 
experimental results is due to the absence of high-resistance 
electrode layers in our theory, which decreases the MR ratio. 
In the experiment, the total RA of the electrode layers is 
estimated to be about lOOmfJ/im^. If we take into account 
this extra resistance, the MR ratio goes down to 10% with the 
RA of a few hundred mil firn^, which quantitatively agrees 
with the one obtained in the experiment. 

The value of the contact radius that maximize the MR ratio 
can be qualitatively explained by considering the effective 
circuit shown in FiglT](b). For the CCP-CPP-GMR spin- valve, 
the resistance of the NOL layer increases with decreasing the 
radius of the contact. In the effective circuit model shown in 
FigUI (b), the change of the contact radius can be realized 
by the corresponding change of effective resistance. Since the 
thickness of the N-layer ijv is much smaller than the spin 
diffusion length Xn, we can employ the following approxi- 
mations: rM{P) ^ f]^ — ppft]^ /{2(P) and rN{AP) ^ oo, 
where d is the contact radius. Based on this consideration, 
in Fig. |4] (b), we plot the MR ratio against the ratio of the 
effective resistances rp/rN- It is observed that the MR ratio 
is maximized when the effective resistances rp and fjv are 
balanced with each other, i.e., rp ~ f^. 

Finally we shall discuss the effect of the spin-dependent in- 
terfacial scattering which is represented by the spin-dependent 
interfacial resistance in the macroscopic model of the 
CPP-GMR |8|. The boundary conditions for electro-chemical 
potentials at the interface located at z — zq are given by 

^J.siz = zo + 0) - ^isiz = zq-O) = rsfs{z = zo), (16) 

where j| is the z— component of the current density and the 
spin-dependent interfacial resistance are defined as 

^f' = (l+7KV2, r^-' = (1 - 7^-72, (17) 

where 7 is the interfacial spin asymmetry coefficient and rf, 
is the total interfacial resistance. 

In our calculation, the spin-dependent interfacial scattering 
is realized by introducing fictitious layers with thickness ^jr, 
resistivity piR = rf,/tiR, spin polarization of the conductivity 
/3iR — 7, and with the infinite spin diffusion length Air 00. 



We carried out the same calculations as shown in Figs |3] 
taking into account of the spin-dependent interfacial scattering 
at both sides of the contact. The parameters are taken to 
be piR = 200f2nm and 7 = 0.62 IS, M\- The results are 
qualitatively same as those shown in in Figs |3] However, the 
contact radius that maximize the MR ratio decreases due to 
the spin-dependent interfacial scattering. For pcu = 17f2nm, 
for example, the contact radius that maximize the MR ratio 
decreases as 3.0nm 2.3nm. 

IV. Conclusion 

In conclusion, we theoretically study the magnetoresistance 
of a CCP-CPP-GMR system by numerically solving the diffu- 
sion equations for spin-dependent electro-chemical potentials 
in three-dimensional CCP geometry. We show that the MR 
ratio is not a monotone function of the contact radius but takes 
a maximum value at a certain value of the contact radius. We 
also show that the value of contact radius that maximize the 
MR ratio is qualitatively understood by considering matching 
of the effective resistances. 
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